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Abstract. We study spin tunnelling in molecular magnets as an instance of a mesoscopic phe- 
nomenon, with special emphasis on the molecule Feg. We show that the tunnel splitting between 
various pairs of Zeeman levels in this molecule oscillates as a function of applied magnetic field, 
vanishing completely at special points in the space of magnetic fields, known as diabolical points. 
This phenomena is explained in terms of two approaches, one based on spin-coherent-state path 
integrals, and the other on a generalization of the phase integral (or WKB) method to difference 
equations. Explicit formulas for the diabolical points are obtained for a model Hamiltonian. 
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1. Introduction 

The last five years or so have seen a burst of activity in the physics and chemistry of 
molecular magnets. This activity is fuelled to a great extent by the ability of chemists 
to synthesize organic molecules containing magnetic ions such as such as Mn, Fe, and 
Co, which then form molecular solids, in many cases of high crystallinity and homo- 
geneity. More than 5000 different magnetic molecular clusters are known, and about a 
100 of these show some behaviour or quirk that is of interest from the persepctive of 
physics [1,2]. These molecules are intermediate bewteen simple paramagnetic salts such as 
CuSO4.K2SO4.6H2O orCe 2 Mg3(N03)i2.24H 2 0, in which the magnetic entites are single 
transition metal or rare earth ions, and superparamagnetic particles of submicrometer size. 
Thus, they show phenomena such as hyeteresis at the molecular level [3-5], but since the 
magnetic entities are molecules and thus well defined with identical size, shape, and orien- 
tation, and since the interactions of the magnetic and nonmagnetic degrees of freedom are 
relatively well understood, there is the potential for understanding the energy relaxation 
process in great detail. Some of these issues are similar to those that arise in the study 
of magnetization reversal of small magnetic particles, and it is hoped that the molecular 
systems will offer insights into the latter, which has obvious implications for magnetic 
recording and storage technologies. 

Our purpose in this article is rather different. The molecular systems have total spin of 
the order of 10, and magnetocrystalline anisotropics of few tens of Kelvin in energy. The 
quantum mechanical dynamics of the spin are then profitably viewed in terms of tunnelling 
in many cases, and it is of interest to look for such behavior. At present, the greatest interest 



1 



Anupam Garg 



is in the molecule [Fes02(OH)i2(tacn)6] 8+ (henceforth abbreviated to Fe§), which is a 
molecular complex in which antiferromagnetic interactions between the Fe 3+ ions within 
a molecule lead to a ground state with a total spin of 10. This molecule shows clear 
evidence for spin tunnelling in the form of a tunnel splitting that oscillates as a function 
of a static external magnetic field [6]. One way to understand this effect is in terms of 
the spin-coherent-state path inetgral for spin. The kinetic term in this path integral has 
the properties of a Berry phase, which can allow for interference of spin trajectories [7,8]. 
In fact, using this property, the oscillations were found theoretically [9], without knowing 
of the existence of Feg. However, another way to understand the effect, which has the 
advantage of using only classical methods of analysis, is in terms of a discrete analogue of 
the phase integral or WKB method. 

In Sec. 2 we shall give a brief discussion of the experiments on Fe§. and in Sees. 3 and 
4, we shall describe the instanton and discrete phase integral (DPI) methods for calculating 
spin tunnel splittings. 

2. Summary of experimental results on Fe 8 

Electron paramagnetic resonance and other measurements show that at low temperature, 
the Fe 3+ spins in one molecule of Fe 8 behave as a single large spin of magnitude 10, with 
an anisotropy energy that is well described by the Hamiltonian, 

Ho = fci Jl + k 2 Jl (1) 

with J = 10, fci w 0.33 K, and k 2 w 0.22 K. The g-factor of the net spin is very close to 
2, with very little anisotropy in the g-tensor. 

If we think of Eq. (1) as a classical energy function for a classical vector J, we see that 
the energy is a minimum when J|| ± z, and a maximum when J|| ± x. We refer to x, 
y, and z as the hard, medium, and easy axes respectively. When the problem is treated 
quantum mechanically, the quantum states corresponding to the two degenerate minima at 
±z should be able to mix together via tunnelling. This picture remains true when we add 
a magnetic field in the xy plane. The Hamiltonian is now 

H = -k 2 J 2 z + (fci - k 2 )J 2 x - g^ B {J x H x + JyHy), (2) 

where we have subtracted out a constant k 2 3 ■ J. The minima are now moved off the ±z 
axes toward the equator, but they continue to be degenerate, and should mix by tunnelling. 
Indeed, as H±, the magnitude of the field in the xy plane, increases, both the angle through 
which the spin must tunnel and the energy barrier decrease, and we expect that the splitting 
A will increase. 

At this point, two questions arise. First, how should one think of the tunnelling of a spin? 
And second, how large is the splitting? One way to answer the first question is to regard 
the last three terms as perturbations that give rise to transitions between various Zeeman 
levels or eigenstates of J z . As usual, we denote the J z value by to. The term gives 
rise to Am — 2 transitions, and thus mixes to = —10 with to = +10 via the —8, —6, 
. . ., +8 states. The H x and H y terms give rise to Ato = 1 transitions, and mix to = — 10 
with +10 via all intermediate levels —9 to +9 (see Fig. la). This picture allows us to think 
of spin tunnelling in direct analogy with a particle tunnelling through an energy barrier. 
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(a) (b) 
Figure 1. (a) Zeeman levels of Feg, showing Am = 1 (solid lines) and Am = 2 

(dashed lines) transitions, (b) The Landau-Zener-Stiickelberg process. 



It is further obvious that if we also apply a field along the z direction so as to tune the 
energies of the —10 and +9 states to resonance, we can also think of tunnelling between 
these states. More generally, we can consider tunnelling between a state with m = m, on 
the negative m side and m — mj on the positive m side. 

The second question is in some sense the main subject of this article. Yet it is useful to 
have an approximate answer before one embarks on a detailed calculation. For a massive 
particle tunnelling in a symmetric double well, the splitting can quite generally be written 
in the form 

/ a \ 1/2 

A = „„„(§) e-\ ,3) 

where u>o is the classical small oscillation frequency about the minima, So is the tunnelling 
action or WKB exponent, and c\ is a constant of order unity. Further, for a smooth poten- 
tial, So can be written as C2V0/W0, where Vb is the energy barrier, and c 2 is another constant 
of order unity, generally close to 5. For the quartic potential V{x) = Vo(x 2 — a 2 ) 2 /a 4 , 
e.g., Ci = 4:^/3 and C2 = 16/3. To apply this approximate formula to the spin problem, we 
may take Vq = k 2 J 2 , but we still need ojo- This, however, can be found by noting that the 
classical vector J can be given dynamics via Hamitlon's equation and Poisson brackets: 

f = < J '*>™ = - J *§- < 4 > 

Then, J x = 2k 2 J y J z , and J y — —2k\J x J z . Linearizing around the equilibrium state J = 
Jz, we find that lu = 2{kik 2 ) 1 ^ 2 J . For Fe 8 , lu 5.4K, and if we use the approximate 
formula (3) with the same c\ and C2 as for the quartic double well, we find that A k> 60 nK. 
This is an extraordinarily low splitting, of the order of 1 KHz in frequency units. At first 
sight, the question of how to experimentally detect such a low splitting would appear to 
be entirely moot, since each molecule experiences stray magnetic fields H' of abcond- 
mat.garg.91 12out200 Oe of dipolar and hyperfine origin. Consequently, the Zeeman levels 
m = ±10 are shifted from perfect degeneracy by an amount e ~ g\i\,JH' ~ 0.1 K, which 
is enormous compared to A. Now to see any resonance at all between two states, the 
energy bias between them must be comparable to or less than the assumptive tunnelling 
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Figure 2. Measured splittings [6] for Fe 8 for (a) —10 <-> 10 transitions for various 

orientations of H in the xy plane, and (b) for H||x between the states m = — 10 and 

m = 10 — n. Figure courtesy of Dr. Wernsdorfer. 



amplitude. In Fe 8 , exactly the opposite is true, and it follows that left to itself, a molecule 
will almost never tunnel at all. 

Wernsdorfer and Sessoli circumvent this problem by applying an ac magnetic field along 
z, in the shape of a triangular wave of period r and amplitude H . As H z changes with 
time, the energies of the to = ±10 states move in opposite directions, and cross at some 
point in the cycle. This induces what are known as Landau-Zener-Stiickelberg (LZS) tran- 
sitions [10,11]. In the limit where the crossing point is passed rapidly, the —10 <-» 10 
transition probability 7 in each passage is very low, a A 2 / (dH/ dt). Since the stray fields 
are fluctuating randomly, different passages may be assumed to be independent, with no 
phase correlation between successive passages. The net LZS rate, or transition probability 
per unit time is given by 



Tlzs = -7 

T 



7rA 2 



(g/j,Bh)AmH 



(7 « 1). 



(5) 



Wernsdorfer and Sessoli measure this rate by measuring the rate at which the magneti- 
zation, after having being initially saturated along the — z direction, say, relaxes to its 
equilibrium direction under the influence of the ac longitudinal field, and a dc field, which 
may or may not have a longitudinal component. 

The outcome of the experiments is shown in Fig. 2. (The temperature at which these 
data are taken is 0.36 K.) Our earlier expectation that A will increase monotonically as H± 
increases is seen to be false. When H||y (<f> = 90°), the behavior is monotonic, but when 
H||x {(f> = 0°), one finds that A oscillates with H x \ 



3. Instanton approach to spin tunnelling 

3.1 General comments on instantons 



For massive particles, the instanton approach is an elegant and effective way to calculate 
tunnel splittings, especially between ground states [12,13]. Extended to spin, the basic idea 
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is to examine the imaginary time transition amplitude 



U 2 i = (n 2 |e 



-HT 




(6) 



where ni j2 are the minima of the classical energy 



E cl (n) = (n|W|n) 



(7) 



and |n) is a maximal-projection spin coherent state, i.e., an eigen state of J • n with eigen- 
value J. 

In the T — > 00 limit, only the lowest two energy eigenstates contribute in a spectral 
decomposition of Eq. (6), and since these differ in energy by A, we get 



ignoring inessential prefactors. If we can calculate U21 in this limit, comparison with 
Eq. (8) will yield A. This calculation is done by appealing to the spin-coherent-state path 
integral 



The paths n(r) all run from ni at — T/2 to n 2 at T/2, and 5[n(r)] is the imaginary time 
or Euclidean action for the path. (This is why the exponent is —5 rather than iS/h.) 

Since T is tending to 00, we can evaluate Eq. (9) via the steepest descents approximation. 
The dominant paths, known as instantons, are just those that minimize the action, i.e. they 
are solutions to the classical equations of motion. The simplest such paths consist of a 
single transit from fii to n 2 . If the scale over which E c \ varies is V, then it follows from 
Eq. (4) that the time scale for this transit is to <~ J/V. For the Hamiltonian (1), e.g., this 
time scale is uJq 1 <~ J _1 (fcifc 2 ) -1 / 2 <C T. Hence, the spin spends most of its time near 
the end points fii , 2 , and the actual transit takes place in a very short time interval. (Hence 
the name instanton.) In the same way, we can find anti-instantons, solutions that go from 
n 2 to ni. We now note that Eq. (4) is autonomous, i.e., does not depend on r explicitly. 
Therefore a translation of the center of the instanton yields an equally good classical path. 
Secondly, one can fit an arbitrary number of instantons followed by antiinstantons into the 
interval T. One finds that the n-instanton contribution to U21 is proportional to T™, and 
the full series is that of a sinh [13]. The A which is obtained in this way can be written as 



where 5; ns t is the action for a single instanton path, and D is a prefactor arising from doing 
the path-integral over small fluctuations about the instanton trajectory. If more than one 
instanton exists, we must add together the corresponding contributions from all of them. 

In the application of this formalism to any tunnelling problem, it turns out that the cal- 
culation of Sinst is often very easy, while the calculation of D is quite hard. (Indeed, 
it turns out that the spin analogue of the van Vleck determinant for the general massive 
particle propagator has only recently been found in a completely logically consistent way 
[14-17].) Since the most interesting aspects of the Feg problem can be understood at the 
level of the instanton action alone, and since the absolute scale of A can be found from the 
DPI method, we shall forego the calculation of D in this article. 



U 2 i ~ sinh(AT), 



(8) 




(9) 



A = Dcxp(~S inst ) 



(10) 
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Figure 3. The kinetic term in the spin action, and its variation. 



3.2 Action for spin 



The interesting features of the spin problem lie in the form of the action, 



/_ 



T/2 



S[n(r)] = iJA[n{r)] + 



E cl [h( T )]dT. 



(11) 



-T/2 



The term _4[n(r)] is the kinetic term, and has the mathematical structure of a Berry phase. 
We write it as 



by which we mean that ^4[n(r)] is the area of the patch S on the unit sphere, whose 
boundary C is the closed curve formed by the path n(r) and a reference path n# (taken 
backwards) running from ni to fi2 (Fig. 3). The reference path ha is arbitrary, but must 
be the same for all n(r) in the path integral. Its choice is equivalent to fixing the gauge. 
However, since we have defined A as a geometrical quantity, an area, it does not depend on 
how we choose coordinates on the unit sphere, and it is obviously nonsingular. Note also 
that on the sphere, each closed curve divides the sphere into two areas, and one cannot say 
which is inside or outside. (This is easily seen by thinking of the seam on a tennis ball.) 
The two areas differ by 4*7r, however, so e~ tJA is identical for the two of them for either 
integral or half-integral J, and hence this ambiguity does not affect the path integral. 

Instead of deriving Eq. (1 1), we shall show that it is correct by checking that its variation 
leads to the classical equation of motion. Suppose we vary the path n(r) to n(r) + SH(t). 
5A is given by the area of the thin sliver enclosed between the curves n(r) and n(r) + 
<5n(r) (Fig. 3). The part of this area due to the segments between r and r + St equals 





(12) 



A(5A) = 5h(r) x [n(r + At) - n(r)] • n(r) 




(13) 



Adding up the contributions from all the segments, we find the total change 
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6A = J Sh(r) ■ x n(r)j dr. (14) 

The variation of the second term in Eq. (1 1) is trivial, and SS can be written as an integral 
of the form J Sh(r) • X where X depends on n and E c \. The extremal condition SS = 
is thus 

ijf x n(r) + ^ = 0. (15) 
dr on 

Taking the cross product with n, and using the fact that n • (dh/dr) = 0, we get 

dn / dE cl \ 
iJ-r- = - I n x — - . (16) 



dr \ dh 

This is exactly what we would get from Eq. (4) with J = Jh and the Wick rotation 
t — > — it. In other words, it is the imaginary time equation for Larmor precession in the 
effective magnetic field dE c \/dh. 

The actual evaluation of A is often more easily done by using Stokes's theorem to trans- 
form Eq. (12) to a line integral. Using notation borrowed from electromagnetism, let us 
write d£L = B • nds, where B(n) = n, and ds is an area element. The line integral is 
§ c A • dh, with B = V x A. Since B = n, it is a monopole field, which, as is known, 
cannot be represented in terms of a nonsingular vector potential. If this singularity is con- 
centrated into a Dirac string at the south pole, we can write 

A[n{r)] = <j> [1 - cos d(4>)] dtf>, (17) 

where C is viewed as parametrized by <f>. This formula is correct as long as C does not pass 
through the south pole, and provided one increments or decrements <j> by 2w every time 
one crosses the date line. 



3.3 Application to Fes 

We now apply the above formalism to a discussion of the ground state tunnel splitting 
in Fe§ when H||x, further limiting ourselves to finding the action, as this is sufficient to 
determine the field points where the splitting is quenched. 

The main problem is to calculate the action for an instanton. To do this, however, we do 
not need to find the actual time dependence n(r), since the energy is conserved along an 
instanton. To see this, note that by Eq. (16) 

dE c \ _ dE ci dh 
dr dh dr 

i dE cl / dE cl \ 
= J^h-{ nX ^h) 

= 0. (18) 

Thus, using energy conservation, we can find the instanton orbit, i.e., the curve traced out 
on the unit sphere #(</>), without finding the actual time dependence. Using Eq. (17), this 
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Figure 4. Interfering instanton trajectories, and numerically computed — 10 <-» +10 
tunnel splitting for a model for Feg with H||x. 



is enough to find the action, as we can take E c \ itself to be zero by adjusting the zero of 
energy appropriately. One point to note is that since fii and fi2 are minima of E c \, there can 
be no real curve connecting these points on which E c \ is the same. The only way to find a 
solution to the instanton equations of motion is to let n become complex. Correspondingly, 
the area A must also be defined on the complexified unit sphere. 

The above discussion is applicable to any tunnelling problem. What is special about Fes 
is the existence of two interfering instantons. If we choose the polar axis to be x (not z), 
and measure the azimuthal angle in the yz plane from y, then we can write the energy as 

E cl (6,<f>) = hJ 2 (cos 9 - cos 6> ) 2 + k 2 J 2 sin 2 9 sin 2 </>. (19) 

We have defined cos# = H/H c , with H c = 2kiJ/gp,B, and added a constant to E c \ so 
that E c \ — along the instanton as discussed above. Writing cos# = u a, the solution of 
this equation gives 

cos = 2 . (20) 

1 — A sin 

However, it is clear that in this equation, we may take <fi to lie in either of the two intervals 
(0, 7r) and (0, — n). Indeed, it is obvious from symmetry that there are two instanton paths, 
which wind about x in opposite directions (see Fig. 4), and Eq. (20) describes both of them. 
If we denote the two paths by A and B, then the imaginary parts of their actions (= i J A) 
are necessarily unequal, since by the interpretation of A as an area, 

S B - S A = i J x fl, (21) 

where f2 [see Eq. (12)] is the area enclosed between A and B (the peanut shape in Fig. 4). 
Using the line integral form (17), we obtain 

0= f (l ^L_W 27 r(l-^=y (22) 

J-n V 1 - Asin 2 0y V \ \/T^V 

The real parts, on the other hand are equal and given by 
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S R = ReS A 



/■±7T 

B = JX 1/2 / 

Jo 



sin 0(1 — Uq — A sin 2 



Nl/2 



J 



In 



v/l-itg + VA" 



1 — A sin 2 

u ,_ ( Vi 1 ~ "oX 1 ~ X ) + u oVX 



:ln 



VAi \7T^A I v /(l- U 2 )(l-A)-w \/A / 



(23) 



The prefactors D are equal for both paths by symmetry, and so, adding together the two 
contributions, we get 



A = 2Dcxp(-S R )cos(Jfl/2), 



which vanishes whenever 



H 



J 



(24) 



(25) 



where n = 0, 1, . . . , 2 J — 1. 

In Fig. 4 we also show the results of an explicit numerical diagonalization of the 21 x 21 
Hamiltonian matrix for our model Hamiltonian for Feg showing clearly that the effect is 
genuine. It should also be clear that the quenching phenomenon is a general one that will 
occur as long as H has the requisite symmetry, and that the detailed form will only change 
the locations of the quenching points. The minima in the tunnelling rate seen by Werns- 
dorfer and Sessoli are spaced approximately 50% further apart than implied by Eq. (25). 
These differences are well understood in terms of higher order anisotropy corrections to 
the model Hamiltonian (1). These corrections do not change the centers of gravity of the 
levels by very much, but they affect the splittings significantly. 



4. Discrete phase integral approach to spin tunnelling 

It will not have escaped the reader that the Wernsdorfer and Sessoli data (Fig. 2b) also 
show an oscillation in the tunnelling amplitude for the m = —10 <-» +9 and —10 <-» +8 
transitions. These oscillations were not predicted in Ref. [9], and are nontrivial instances of 
a diabolical point, or conical intersection [18,19], a degeneracy obtained at isolated points 
in a two dimensional parameter space. In Fe§, the parameters on which the Hamiltonian 
depends are H x and H z . 

In retrospect, the new oscillations can also be understood in terms of path integrals, 
but it is simpler to adopt a different point of view, namely, the discrete phase integral (or 
WKB) method. The basic idea of this method is to directly solve Schrodinger's equation in 
the J z basis as a recursion relation or difference equation, by making use of the similarity 
between difference and differential equations, and exploiting WKB type ideas used to solve 
the latter. 

To see what is meant, suppose that \ip) is an eigenstate of H with energy E. With 

J z \m) — m\m), {m\ip) — C m , (m\H\m) — w m , and (m\H\m') = t m . m < (m ^ m'), we 
have 

w m C m = EC mi (26) 

n 
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where the prime on the sum indicates that the term n — m is to be omitted. This equa- 
tion can be thought of as a tight binding model for an electron in a one-dimensional lat- 
tice with sites labelled by m, and slowly varying on-site energies (w m ), nearest-neighbor 
(i TOjTO ±i) hopping terms, next-nearest-neighbor (i TO>m ±2) hopping terms, and so on. Un- 
der certain conditions (which we shall see below) we can understand the motion of the 
electrons in terms of wavepackets which are well localized in both quasi-momentum and 
real space. This is of course the approximation of semiclassical electron dynamics. As 
discussed in standard texts on condensed matter physics, this approximations works when 
the wavepacket can be assigned a real space width Am much less than the length scale 
over which the properties of the electron band (in our case the w m and t m ^ m+a ) vary, and 
at the same time, a momentum space width Aq much less than the reciprocal bandwidth 
2tt. Since the two widths must be constrained by the uncertainty principle AqAm <~ 1, 
the real space width Am must be much greater than the lattice constant (in our case 1). 
Formally, we need to be able to find smooth continuum approximants w(m) and t a (m), 
such that 

w(m)=w m , (27) 
t a (m) = (i TO , m+a + i m , m _ a )/2, a = 1,2,... (28) 

whenever m is an eigenvalue of J z , and that these quantities vary slowly enough that 

dw_ = / w(m) \ dta_ = / t a (m) \ ^ 
dm \ J J ' dm \ J J 

These conditions are generally met when J> 1, which is also the semiclassical limit for 
spin that we expect intuitively. 

Given the conditions (29), the basic approximation, which readers will recognize from 
the continuum case, is to write the wavefunction as a linear combination of the quasiclas- 
sical forms 

(i I" q(m')dm'), (30) 



where q(m) and v(m) obey the equations 

E = w{m) + 2 ^ t a (m) cos(a(7) = H sc (q, m), (31) 

a>l 

v(m) - dH sc /dq. (32) 

Equations (31) and (32) are the lattice analogues of the eikonal and transport equations. 
Equation (30) represents the first two terms in an expansion of log C m in powers of 1/ J. 

Previous work with the DPI method [20-24] has been limited to the case where the re- 
cursion relation has only three terms, i.e., only nearest neighbor hopping is present. Braun 
[23] mentions many problems in many areas (nuclear and atomic physics, to name just 
two) where the Schrodinger equation turns into a three-term recursion relation in a suitable 
basis. All the types of problems as in the continuum case can then be treated — Bohr- 
Sommerfeld quantization, barrier penetration, tunneling in symmetric double wells, etc. 
In addition, one can also use the method to give asymptotic solutions for various recur- 
sion relations of mathematical physics, such as those for the Mathieu equation, Hermite 
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polynomials, Bessel functions, and so on. The interesting features all arise from a single 
fact — that the DPI approximation breaks down at the so-called turning points. These are 
points where v(m) vanishes. One must relate the DPI solutions on opposite sides of the 
turning point by connection formulas, and the solution of all the diverse types of problems 
mentioned above depends on judicious use of these formulas. 

For the Fes Hamiltonian (1) with arbitrary H, the recursion relation involves five terms. 
The diagonal terms (w m ) arise from the and J Z H Z parts of H, the t TO , m ±i terms from 
the J X H X and J y H y parts, and the t m ^ m ± 2 terms from the part. This seemingly minor 
modification is responsible for all the beautiful spectral features seen experimentally. 

In the three term case, only w(m) and ti(m) are present, so H sc (q,m) = w(m) + 
2ti(m)cosq, and v(m) = —2smq(m)ti(m). Hence, the turning points arise when q = 
or q = ir, corresponding to the local m-dependent band edges. This is completely 
analogous to the condition V(x) = E that defines a turning point for a massive particle 
moving in a potential V(x), in that a turning point is a limit of the classically allowed range 
of motion. 

In the five term case, by contrast, 



Hsdq, to) = w(m) + 2t\{m) coaq + 2t2(m) cos(2g), (33) 



so that 



v(m) = — 2smq(m)(ti(m) + 4t2(m) cos q(m)). (34) 
In addition to q = and q = n, now the velocity can also vanish when q = q*(m), where 

cosq»(m) = -ii(ra)/4t 2 (m). (35) 
If we make use of the eikonal equation, we see that a turning point arises whenever 

E = U (to) , U* (m) , or U„ (m) , (36) 

where, 

U (m) = H sc (0, to) = w(m) + 2t 1 {m) + 2t 2 (m), (37) 

U n (m) = H sc (7t,to) = w{m) - 2t 1 (m) + 2t 2 {m), (38) 

U,(m) = H sc ( qif ,m) = w(m) - 2t 2 (m) - t .^\ . (39) 

4t2(TO) 

We call these curves critical curves; collectively they play the same role as the potential 
energy V(x) in the continuum phase integral method in determinimg the turning points. 

The new turning point (m) and the corresponding critical curve {/* (to) have no ana- 
logue in the continuum case. This turning point need not be at the limit of the classically 
allowed energies for fixed to, and may lie strictly inside the band, or in the forbidden re- 
gion outside it. We show the general form of the critical curves for Fe 8 when H has both 
x and z components in Fig. 5. The point m' c is of this new type, and here it may be said to 
lie "under the barrier." 

From this point on, the analysis of eigenfunctions and eigenvalues is conceptually 
straightforward, although quite lengthy, and we refer readers to already existing papers 
for the details [25,26]. The basic idea is simple. For an energy E such as that shown in 
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Figure 5. Critical energy curves for the Feg Hamiltonian when H has both x and 
z components. In the right hand figure, the region near the left minimum of Uo is 
magnified, showing the various turning points. 



Fig. 5b, the region m' a < in < m' b is classically allowed, along with its counterpart for the 
right hand well. To find the eigenfunctions, one starts with a solution that is exponentially 
decaying to the left for m < m' a , continues it through the turning points m' a , m' b , and m' c 
by the use of connection formulas at each point, and matches it with a similar solution 
obtained by starting from the right. In general the wavefunctions will agree in the central 
region near m — only if the energy E is properly chosen — this gives the eigenvalue 
condition. Nevertheless, as already stated, the calculation has several novel aspects. For 
example, refering to Fig. 5b, one should clearly have Im q(m) < in the region m < m' a . 
In conventional WKB, there would be one solution satisfying this demand, but now there 
are two. As another example, the two solutions in the region m' b < m < m' c which are 
both decaying exponentially to the right, turn into decaying solutions with an oscillatory 
envelope as the point m' c is crossed. It is this oscillatory exponential solution (along with 
the concommitant fact that there must be two such solutions when they exist) that allows 
the tunnel splitting to vanish if H x and H z are chosen properly. Cutting short an already 
long tale, we give the results for the location of these points. The point where the £'i\\ level 
in the negative J z well (with £' = being the lowest level) and the £"th level in the positive 
one are degenerate, is at H y = 0, and 

H c 2J y ' 

S-&n = £p [J - a - W + t . + 1)h (4i) 

with n = 0,1, ... ,2J -(£'+£" + 1). Here, A = k 2 /k lt and H c = 2k 1 J/g^i B . Exactly 
as seen in the experiments, the — 10 <-> 9 points are shifted by half a period with respect to 
the — 10 <-» 10 points, the — 10 <-> 8 points are shifted by another half-period, and so on. 
Further, many of the degeneracies occur simultaneously, i.e., at the same values of H x and 
H z . This latter feature has only been tested indirectly in the experiments so far. 

We conclude with an amusing fact about Eqs. (40) and (41). These are the leading 
1/ J results from a semiclassical analysis. Nevertheless, for the model (1), they are exact 
as written! This has been demonstrated in Ref. [27], and while the exactitude is spoilt by 
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the higher order anisotropics present in real Fe§, the exact results may pave the way for a 
quantitative treatment of this and other perturbations. 
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